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Abstract 

Within the framework of the minimum quadratic Poincare gauge theory 
of gravity in the Riemann-Cartan spacetime the dynamics of homogeneous 
anisotropic Bianchi types I-IX spinning-fluid cosmological models is investi- 
gated. A basic equation set for these models is obtained and analyzed. In 
particular, exact solutions for the Bianchi type-I spinning- fluid and Bianchi 

type-V perfect-fluid models are found in analytic form. 
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Typeset using REVTpX. 



I. INTRODUCTION 



As is well known in the modern theories of fundamental physical interactions the princi- 
ple of local gauge symmetry takes the most important position. On the basis of this principle 
some unified gauge theories - the Weinberg-Salam theory, quantum chromodynamics, Grand 
Unified theories, etc. - that successfully describe electromagnetic, weak, and strong inter- 
actions of elementary particles, are constructed and carefully studied. Therefore, it seems 
natural and attractive to construct the theory of gravity as the gauge theory of spacetime 
symmetries and to study its consequences. During the last 35 years various versions of gauge 
theory of gravity have been suggested and investigated. These theories differ one from an- 
other by the choice of the gauge group (Lorentz, Poincare, de Sitter, metric-afline, etc.) and 
gravitational Lagrangian (linear, quadratic ...). However, the importance of the Poincare 
symmetry in particle physics lead one to consider the Poincare gauge theory (PGT) (see 
[1-5]) as the most suitable for description of the gravitational interaction. 

In the PGT the gravity is described by two set of gauge fields: the tetrad fields h\ , 
and the Lorentz connection coefficients A lk u , the spacetime continuum is represented by a 
four-dimensional differentiable Riemann-Cartan manifold U4. 

The corresponding gauge field strength tensors are the torsion: 

and the curvature: 

F%, = 2d [M A%+2A%A fc b>] . 

The sources of the gravitational field in this theory are the canonical energy-momentum 
tensor (EMT) V, and the tetrad spin current (TSC) J ifc ". 

Investigation of gravitating systems in the approximation of homogeneous isotropic 
spaces shows, that by including in the gravitational Lagrangian the terms quadratic in the 
curvature and torsion tensors the PGT allows to prevent the appearance of the gravitational 
singularity with infinite energy density [3, 5-8]. In connection with avoiding the cosmo logical 
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singularity problem, the study of more general homogeneous anisotropic models in the PGT 
is of certain interest. 

In this article we shall study the dynamics of homogeneous anisotropic models within the 
framework of the so-called minimum quadratic (Poincare) gauge theory of gravity (MQGT) 
- the simple version of the PGT with the following gravitational Lagrangian: 

L 9 = h(- X + / F + aF 2 ), F F = FV, F = F%, (1) 

where h = deth\, F is the scalar curvature of the Riemann-Cartan spacetime U4, 
fo = (167rG) _1 , (G is the Newton's gravitational constant), x is a cosmological term, 
and a < is a dimensionless coefficient (in the sequel we use all the notations of [9-12] 
unless otherwise stated). As was shown in [5, 9] the MQGT satisfies all restrictions following 
from the simultaneous consideration of the quantization problem (the theory without ghosts 
and tachyons) [13], Birkhoff's theorem [14], and the problem of cancelling the metric singu- 
larity in the homogeneous isotropic cosmological models [6]. Satisfying the correspondence 
principle with General Relativity (GR) for the systems with rather small energy density, 
the given theory gives a chance to overcome some difficulties of GR, which appear on the 
description of gravitating systems under extremal conditions (extremely high energy density, 
spinning matter...). 

The paper is organized as follows. In section 2, for spinning-fluid systems we obtain 
the gravitational equations of MQGT in the synchronous frame of references. In section 
3, we derive the equations of motion for Bianchi models in the MQGT. In section 4 and 5 
we find analytic solutions of the obtained basic equation set in the simple cases: Bianchi 
spinning-fluid type-I and Bianchi perfect-fluid type-V models, respectively. Some interesting 
properties of the found solutions are analyzed in section 6. 
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II. MQGT IN THE SYNCHRONOUS FRAME OF REFERENCE. 



By independent variation of Lagrangian (1) with respect to gauge fields h\, and A %k v 
we obtain two field equations of MQGT in the form: 



2(/ + 2aF)F"i - (/ + aF)Fh" i + = ~t M *, 



2V„[(/ + 2aF)h [i ' t hfc]1 



— T v 



(2) 
(3) 



It follows from (2) that F = ^(t£ + 4%). Then, it is easily to solve eq. (3): 

s V = (i - - M- 1 ( - v A + ^J,^) + ^jjAtf) • 



(4) 



By virtue of (4), after some simple transformations (see Ref. [9]) we can represent (2) 
in the form of the Einstein's equation with the effective EMT T A jj constructed from EMT 
tf and TSC as: 

G A ^({}) = —Tiff, (5) 
47b 

where G A jy({}) is the ordinary Einstein's tensor and 



t a ; 7 = (i - 4(3 X - Kr 1 [t a ^ - + 4 X ) 2 + X g x " + ^f o ( 1 ~ ^ 



+ 



+2/ [v M 9 A M - -d^d x M - g^V^M + -^M^M) 
2 4 



(6) 



where M = ln|l - A(3x - P = ~ 7* > °, tV = tV + (V„ + 2S„) (J' 4 *" + J A ^ + J^ A ), 

Jo 

St- = S a au, and Vj, denote covariant derivatives calculated by means of the total 

connection r A MT and the Christoffel symbols 7 A MT respectively. 

As the source of gravitational fields we consider the perfect fluid with angular spin 
momentum (spinning fluid). In [15-18] a variational formalism for relativistic dynamics of 
spinning fluid was developed in the Riemann-Cartan spacetime. In the simplest spinning 
fluid model the canonical EMT t A ^ and TSC J MZ / take the following form: 



t A ^ = (p + pjuV -pg V +j A r u M u a V CT u r , 



(7) 



where p is the invariant energy density, p ia the pressure, u M is the 4-velocity vector, 
= — is the spin density tensor. By using the relations: ]^ v n u = from (7) we get 
t£ = p — 3p . Then, the torsion tensor (4) can be rewritten as 

SV = (1 " ~ »- : ( - ^u A + ^j^tf). (9) 

In view of (7)-(9) and the equations of rotational motion (see [15,16]) the effective EMT 
T^j can be transformed to the form: 

* 2 

= (l-4/3 X -^)- 1 [(p + p-^-(l-4/? X -^)- 1 )u A u^-g A ^p- X + 

+ A X f - i(l - 4/3 X - /ft*)" 1 ) + j"<V>ft,M - u"u T u("V„j A > T - 
-V^V*)] + 2/ [v^<9 A M - i(9 M M(9 A M - g^V^M + \dMd v M)] , (10) 
where j 2 = \}^f u . 

Following [19] we rewrite the MQGT gravitational equations in the synchronous frame 
of reference, in which the spacetime metric takes the form: 

ds 2 = dt 2 - 7Q/3 dx Q dx /3 , (11) 

where t is a synchronous time, x a — spatial coordinates (first greek letters a, (3 ... pass 
values 1, 2, 3). As in [19] let us introduce the notation: 

7~ta.fi = la(3 (12) 

(a dot denotes differentiation with respect to the synchronous time t). The quantities Hap 
make up a 3-dimensional tensor, raising or lowering of one's indices, and so the covariant 
differentiation is calculated by the 3- metric 7 Q( g. In terms of 7^,3 and H. a p the Einstein 
tensor's components can be expressed as: 



1 1 * 1 



0" = \^aHi-VpHi), 



Gl = \ni + \inZ + vl - l -5i{nl + \ini + \n]n s , + vf), (is) 

where 7 = det^ap), is a 3-dimensional Ricci tensor constructed from 7 Q/3 . Concerning 
the source of gravitational field we will consider only the case of rest matter, when u° = 
1, u a = 0, (hence we have j M0 = 0). Using these relationships we obtain the following 
expressions for the components of T^ff i n the synchronous frame of reference: 



0(e//) = (l-A(3x-f3t:)-'[p + x-jK + AxY-—i 2 (l-Af3x-f3t:)- L ]- 

-2/ (^MK-^M 2 ), 



L <*(eff) 



1 a(e//) 



-(1-4/3* -/ft") 



a\-l 



--(l-APX-M-'Vpfa, 

[p-x+jK+ 4x) 2 - ^j 2 (i - wx - ^r'n- 



\inif, + n^\)] + Vo[\mK + ( - m + Imw} + \u 2 )5i 



(14) 



Note that covariant differentiation V of the quantity in the expression of ^a(eff) 
is calculated by the 3- metric jap. 



ILL BIANCHI MODELS IN THE MQGT: A BASIC EQUATION SET 

Let us proceed to derive the equations of motion for homogeneous anisotropic (Bianchi) 
cosmological models. Homogeneity of the 3-space implies invariance of the spatial metric: 

dl 2 = 7a/3 dx a dx /3 

at each fixed moment of time t with respect to the 3-dimensional group of motion G 3 
(see [19-21]). So that local characters of homogeneous anisotropic models are determined 
fully by types of the corresponding motion's group G, which was classified by Bianchi with 
respect to independent sets of the structure contants C c a b . Here and later on, the first latin 
letters a, b,c... take the values 1,2,3 and number of the basis's vectors. 
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Let X a = e a a ^, (a = 1, 2, 3) be generators of G, that satisfy the following commutation 
relations: 

[ X a' X fc] = C C afe X c , (15) 

where e a a are the basis's vectors of a 3-dimentional space with the metric ^ a p. As we see 
from (15) the structure constants are antisymmetric in their lower indices: G c ab = — C c ba , 
and satisfy the Jacobi identity: 

Ce /-id I (~\e r~\d , r\e /~\d n 
ab^ ec -r \~> bc^ ea "T O ca^eft — U. 

We also introduce basis's vectors e£ that satisfy the relationships: 

« = #5 « = 51 (16) 

Expanding the spatial metric with respect to the basis (e a a ) we obtain the following 
time-dependent functions: 

Vab(t) = la^a^b- (17) 

The Bianchi classification of homogeneous spaces reduces to determining all nonequiva- 
lent sets of structure constants. In place of the three-index constants we introduce a set of 
two-index quantities, defined by the dual transformation: C c ab = e abd C dc , ( e abc = e abc is 
the unit antisymmetric symbol with £123 = 1. Following [22], we decompose nonsymmetric 
C ab into a symmetric tensor n ab and a vector a c as: 

Q ab = n ab + £ abc ac 

It is easy to diagonalize the symmetric tensor n ab ; let n 1; n 2 , n 3 be its eigenvalues. It 
then follows from the Jacobi identity that if a a exits it is an eigenvector of n ab corresponding 
to a zero eigenvalue: 

n afe a a = 0. 

Without loss of generality the vector a a can be written CIS £1(2 (^J 

0, 0) with nia = 0. 

Then, we get either n x or a is equal to zero. Bianchi models with a = were defined by 
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Ellis and MacCallum as A-class models, also known as tensor models. Bianchi models for 
which n x is set equal to zero are known as B-class models - vector models. By rescaling 
the operators Xi, X 2 , X 3 one can set the eigenvalues ni, n 2 , n 3 to be equal to 0, ±1. The 
classification of these models is given in the tables A, and B, respectively (see [19-20]): 



Table A: Classification of A-class models (a = 0) 



Type 


ni 


n 2 


n 3 


I 











II 


1 








VII 


1 


1 





VI 


1 


-1 





IX 


1 


1 


1 


VIII 


1 


1 


-1 



Table B: Classification of B-class models (n x = 0) 



Type 


a 


n 2 


n 3 


V 


1 








IV 


1 





1 


III 


1 


-1 


1 


VI h 


a > 


-1 


1 


vn h 


a > 


1 


1 



By virtue of (12) and (15)-(17) the expansions of H%, Pf, (V 'j£ - VpWg), and V#j£ 



with respect to {e a a } lead to the following expressions (see [19]): 

HI = Hie\e\ = rj ac r, bc , H a a = H%, 



1 

277 



2C M C ad + C db C ad + C hd C da - C d d (C b a + C a 6 ) + S b a [(C d d ) 2 - 2C cd C cd ] 
e Q a(V a W? - VpHi) = n c b (C b ca -5 b a C d dc ), 



(18) 

(19) 
(20) 



e a a(V j aj /3 a ) — -fcC a b C — j b a C C bc, 
where r\ = det(r] ab ), ] h c = j /3 7 e 7 c e fe / g. and, in accordance with the general rule: 



(21) 



C° a = Vac C cb ; C ab = VacVbd C 



cd 



Substituting (17)-(21) to the expressions (13), (14) we find the following expansions of 
the components of G A ^ and with respect to the basis {e a Q } : 



GO uc((~\b Kbi~\d \ 



4 ?7 



(22) 



rpO 

1 0(e//) 



i a(e//) 



~0- - 4/3x - ^)" 1 (d 6 c C ab c - j 6 a C c fec ), 



rpb 

L a(eff) 



- (1 - 4/3 x - pt^r 1 [p - X + jK + 4 X) 2 ] + 57(1 - 4/3x - K)- 2 
-tC 1 - 4 ^ - « _1 (jc 6 ^ + J c a^) - 2/ [(M + ±MW« + ^M 2 )^ - \m-HI 



(23) 

Thus, we obtain the MQGT basic equation set describing the dynamics of Bianchi (types 
I-IX) spinning-fluid cosmo logical models, written in the basis {e Q a }, in the following form: 



G° = 



1 rpO 

w ; o(e//) ' 



pO _ ^ rpO 

a ~~ 2fe a(e//) ' 



(24) 



where the corresponding compononents of G Am and T A ^ are determined by the expressions 



(22)-(23), which contain only the time-dependent functions. 

Note that these equations reduce to the corresponding equations of the Einstein-Cartan 
theory in the case when (3 = 0, and coincide with the ordinary Einstein's equations of GR 



in the case when (3 = and j 



0. 



In the following sections we will study Bianchi models of types I and V, which are the 
direct anisotropic generalization of the flat and the open isotropic Friedmann-Robertson- 
Walker models respectively 
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IV. BIANCHI TYPE-I SPINNING-FLUID MODELS: ANALYTIC SOLUTIONS 



We now study the dynamics of the simplest case of homogeneous anisotropic - Bianchi 
type-I - models, when all of the structure constants C c ab = (see also [12]). Then, the 
metric of 3-space can be written in the diagonal form: 

Vab = diag(rl(t), 4(t), r 2 3 (t)). (25) 

In this case we note that under the presence of spin momentum j a b 7^ the 3-space will 
be either isotropic or axially symmetric. In fact, for a 7^ b we find from (22) G£ = and 
from (23) 

(there is no sum on indexes a, b here). It follows that 

(^-^)j«6 = (a, b= 1,2,3; a ^b). (26) 

In the spinless matter case (j a 6 = 0) relation (19) is automatically satisfied. If at least two 
of the three components of ] a b are not equal to zero, then (26) implies ^ = ^ = ^, which 
leads to the isotropic 3-space: r 1 = r 2 = r 3 . If only one of the three components of j a t, is 
nonzero, for example j 12 , it follows that the 3-space will be axially symmetric ri = r 2 7^ r 3 
By analogy with GR we introduce new variables: r = (r^^) 1 / 3 , h = (r/r), h^ = 
(rjfc/rfc) (k = 1,2.3). Let us consider the general case, when t° = p — 3p 7^ 0. [] Then, it 
is easy to find 

H b a = 2diag(h u h 2 , h 3 ); rj = r 6 ; V b a = 0. (27) 
By using (27) we can rewrite (22) and (23) as 

1^3 2 9 2 

x The special case of radiation t£ = p — 3p = was carefully studied in [11]. 
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G| = K - 3h + 3hh a - ^ELi h ' " ^ 

T° (e//) = (l-4/3 X -^)- 1 [p + x-f(t^ + 4 X ) i 



rpa 



-(1 - 4/3 X - Z^)" 1 [p - X + + 4x) 2 



+ 



+ i-(l - 4/3 X - (3tlY 2 - 2/o [M + M(3h - h a ) + ^M 2 

0/0 L 4 



(28) 



After some manipulations from (24) and (28) we obtain (for the case of x — see also 
[10]) 

h k = h + Sfe r- 2 (l-4/3 X -^)- 1 ; ELi Sfc = °' ( 29 ) 



h-fa-^x-^)- 1 ^) 



;i-4/3 X -/3t: 



a\-l 



4(p + X-^ + 4 X ) 2 ) + 



+ (l-4/3x-^)- 1 r- 6 I 2 - 



48/ 2 



(30) 



h + 3h 2 = (1-APx-PK) 



a\-l 



4/0 

+f((0" + 5h(0) 



(p-p + 2 X -^ + 4 X ) 2 ) + 



(31) 



where s*; (A; = 1,2,3) are integration constants and I 2 = |l]^ =1 s 2 is a measure of the 
anisotropy. 

In the case considered, the spin conservation law (by using the equation of rotation of 
fluid in [15]) takes the following simple form: 



(j)- + 3hj = 0. 



(32) 



and "the energy- momentum conservation law" V^T^ = yields 



p + 3h(p + p) = 0. 



(33) 



From these conservation laws we obtain 



exp\ — 



dp 



3 J p + p(p) 



jr 3 = J = const. 



(34) 
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Assuming p = p(r), j = j(r) in accordance with (34), it is easy to find the solution of 
(30) in an analytic form: 

t-t = / -i^-dr, (35) 

J*0 $ 2 ' (r) 



where 



$i(r) 



3r i"4^-^ + f (p + p)(l- 3 |) 



3 



2/c 



o 



P. 



*2(r) = ^(l-4^x-» P + X-t(^ + 4 x) 2 + 



4 



A 



r 6: 



and A = 9(l 2 - ^jjL). It follows from (29) that 

r fc = r exp(s k J r~ 3 (l - 4/3 X - /ft*)" 1 ^) . (36) 

Thus, under the given equation of state p = p(p) the analytic solutions (34)-(36) ade- 
quately describe the dynamics of the Bianchi spinning-fluid models. Note that the solution 
(35) makes sense only if <3> 2 (r) > 0. For the case of vanishing cosmological term x = the 
study of these solutions (for more details see Ref [23]) shows that we can construct Bianchi 
type-1 cosmological models, which are regular in the metric and in the torsion, by fulfilling 
the condition A < (i.e. when the spin effect dominates over the effect of anisotropy). For 
the case of x an d the linear equation of state: 

P = ^pP, < 7 < 1, (37) 

a complete qualitative analysis of properties of the obtained solutions is perfomed in [12] by 
using methods of qualitative theory of dynamic systems [24]. 

V. BIANCHI TYPE-V PERFECT-FLUID MODELS: ANALYTIC SOLUTIONS. 

In this section we restrict ourselves to the spinless matter case of the perfect fluid, when 
— 0. As was shown in [25] in the general case of Bianchi models (except for type-I) the 
metrics of anisotropic models can not be diagonalized in the presence of spin momentum. 
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Taking (xi, x 2 , x 3 ) as local coordinates we can write the interval of Bianchi type-V 
models in the diagonal form [25, 26]: 

ds 2 = dt 2 -ri(t)dxl-e 2xi (ri(t)dxl-rl(t)dxl). 



For the models considered we have the following structure constants: C 21 = — C 2 2 = 
C 31 = — C? 3 = I. Then, it is easy to find the following relationships: 
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K = 2diag(h 1 , h 2 , ha); V = r b ; V b a = --5 b a . 



(38) 



By using (38) and the restriction: ] ab = from (22) and (23) we find: 



rpa 



I V 3 h 2 - -h 2 + — 
2 Z^=A 2 r 2 ' 

G° = 2h 1 -(h 2 + h 3 ), 
G| = h a - 3h + 3hh a - \j2l =1 h l ~ \^ + ^ 

T^eff) = (i-^x-M-^p + x-^t^ + ^-e/o^Mh+iM 2 ), 

= -(l-4/5 X -«- 1 [p-x+^(t^ + 4 X ) 2 ] -2/ [M + M(3h-h a ) + lM 2 ]. (39) 



It is clear that equation = — 2X^i(e//) = ^ leads to the relationships: 



h = hi = -(h 2 + h 3 ). 



(40) 



It follows from (40) that r ± — r . By virtue of (40) it is easy to show that the equations 
= — 2^^a(e//) nave ^ ne following integrals of motion: 

h 2 = h + 1(1 - 4/5 X - /ft?)" 1 ; h 3 = h - 1(1 - 4/5 X - «-\ (41) 

and consequently 

r 2 = Yexp(sjr- 3 {1 - 4(3 X - /3t*) -1 cft); r 3 = rexp^^! - 4/3* - Pffidi) , (42) 



where s is an integration constant and has the sense of a measure of the anisotropy. By 



virtue of (40)-(41) the equation G° = ~^ r ^o( e ff) can ^ e transformed to 



(l-APx-K)- 1 
13 



4 



+(l-4/3 x -/3t 



^2 n 



<T\-1 



3r 6 



+ 



1 



(43) 



The obtained equation is a direct generalization of the corresponding equation for open 
homogeneous isotropic models in MQGT (see [6]), that is the case, when the measure of 
the anisotropy s = 0. In the given case the " energy- momemntum conservation law" 
V M T^ 7 = has the same form as eq. (33). In a way similar to that of section 4 we can 
find solutions of eq. (43) in an analytic form: 

$s(r) 



t-tn 



ro $l /2 (r) 



dr, 



(44) 



where 



l-4/3 x -^ + ^(p + p)(l-3^) 



*4(r) = ^f(l-Wx-/K)[p + X- j(t^ + 4 X ) 2 ]+i(l-4/? X -^) 2 + 



3r 6 ' 



VI. DISCUSSION 



The exact solutions (37) and (44), that were found in analytic form, describe dynamics 
of the spinning-fluid Bianchi type-1 and perfect-fluid Bianchi type-V respectively. The cor- 
responding torsion functions have been written in (9). A brief analysis shows that the given 
solutions are, in general, "weakly singular" i.e have following properties: 

- an energy density is finite at any time: p < p cr < oo; 

- a "volume" of "3-space" is positive V = r 3 = (rir 2 r 3 ) > V cr > at any time; 

- one or two of the three metric functions r^ may become zero at a finite moment of time 
when 1 - 4(3 X ~ P(p - 3p) = 

- torsion functions diverge S(t) = oo, when 1 — A(3x ~ P(p — 3p) = 



However, in some cases (depending on the values of the spin density, the anisotropy 
measure, the cosmological constant, the parameter (3 , and the equation of state) we can 
find truly regular solutions with finite energy density, nonzero metric functions and finite 
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torsion functions. Using methods of qualitative theory of dynamical systems we can make 
a detailed analysis of properties of every possible solutions of the MQGT equations for the 
models considered. 
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